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We study generation of magnetic field in primordial plasma of the standard model (SM) particles 
at temperature T > 80 TeV much higher than the electroweak scale. It is assumed that there is 
an excess number of right-handed electrons over left-handed positrons in the plasma. Using the 
Berry-curvature modified kinetic theory to incorporate effect of the Abelian anomaly, we show that 
this chiral-imbalance leads to generation of hyper-magnetic field in the plasma in both the collision 
dominated and the collisionless regimes. It is shown that in the collision dominated regime the 
chiral-vorticity effect can generate finite vorticity in the plasma together with the magnetic field. 
Typical strength of the generated magnetic field is 10 27 Gauss at T ~ 80 TeV with the length scale 
10 °/T whereas the Hubble length scale is 10 13 /T. Further the instability can also generate the 
magnetic field of order 10 31 Gauss at typical length scale 10/T. But there may not be any vorticity 
generation in this regime. We show that the estimated values of magnetic filed are consistent with 
the bounds obtained from present observations. 
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Introduction 


There is a strong possibility that the observed mag¬ 
netic fields in galaxies and in inter-galactic medium could 
be due to some process in the very early Universe. Un¬ 
derstanding the origin and dynamics of the primordial 
magnetic field is one of the most intriguing problems of 
the cosmology (see the recent reviews [1-3]). It should 
be noted here that there still exists a possibility that 
the fields may not be of primordial origin but might 
be created during the gravitational collapse of galaxies 
[4, 5]. In the present work we are interested in the pri¬ 
mordial origin of the magnetic fields. There exist sev¬ 
eral models describing the generation of primordial mag¬ 
netic fields in terms of cosmological defects [6-8]; phase 
transitions[9-12]; inflation[13, 14]; electroweak Abelian 
anomaly[15, 16]; string cosmology [17, 18]; temporary 
electric charge nonconservation[19]; trace anomaly [20] or 
breaking gauge invariance [21], In a recent work[22], it 
was shown that the process like Biermann battery can 
play a role in generating the primordial magnetic field 
just after the recombination era. 

In recent times there has been a considerable inter¬ 
est in studying role of the quantum chiral anomaly in 
generation of the primordial magnetic field [23-25]. In 
Ref. [15] (see also [16]) it was argued that there can be 
more right-handed electrons over left-handed positrons 
due to some process in the early Universe at temper¬ 
atures T very much higher than the electroweak phase 
transition(EWPT) scale (~100GeV). Their number is ef¬ 
fectively conserved at the energy scales much above the 
electroweak phase transitions and this allows one to intro¬ 
duce the chiral chemical potentials At tempera- 
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ture lower than Tr ~ 80 TeV processes related with the 
electron chirality flipping may dominate over the Hubble 
expansion rate and the chiral chemical potentials are not 
defined [26-28]. Further the right handed current is not 
conserved due to the Abelian anomaly in the standard 
model (SM) and it satisfy the following equation: 
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here, y lit/ = d l 'Y 1 ' — is the field tensor associated with 

the hyper-charge gauge field Y 1 ' and y = l/2e MI/pA TpA- 
£y and By respectively denote hyper-electric and hyper- 
magnetic fields. Further, g is associate gauge coupling, and 
yn=- 2 represents hypercharge of the right-electrons. The 
right-hand side of the first and second equality signs are re¬ 
lated with the Chern-Simon number Ncs- 
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The anomaly equation (1) relates change in the right handed 
electron density with the variation of the topological (Chern- 
Simon or helicity) charge of the gauge fields. It has been 
shown in Ref. [29], that CS term contributes in the effective 
standard model Lagrangian of the field V M , by polarization 
effect through non-zero mean pseudo-vector current Jj 5 = 
g 2 y 2 Rl 2 < eR 7 j 7 5 e_B >= -g ,2 y 2 R y e RBj/An 2 and the effective 
Lagrangian for gauge field V M in SM is [29, 30]. 
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If these hyper-magnetic fields survives at the time of EWPT, 
they will give ordinary magnetic fields due to the electroweak 
mixing = cos0 w Y/. Where is the massless U{A)y 
Abelian gauge hypercharge field. It was shown that the chiral 
imbalance in the early Universe can give rise to a magnetic 
field B ~ 10 22 G at temperature T ~ lOOGeV with a typical 
inhomogeneity scale ~ 10 6 /T[15]. In this work the authors 
have studied the Maxwell equations with the Chern-Simon 
term and a kinetic equation consistent with equation (1). It 
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was found that the transverse modes can become unstable and 
give rise to the hypercharge magnetic held [31, 32]. In Ref. [23] 
the authors have used magnetohydrodynamics in the presence 
of chiral asymmetry to study the evolution of magnetic held. 
They have shown that the chiral-magnetic [31, 33, 34] and 
chiral-vorticity effects [35] can play a signihcant role in the 
generation and dynamics of primordial magnetic held. Fur¬ 
ther, it was demonstrated in Ref. [25] that evolution of the 
primordial magnetic held is strongly influenced by the chi¬ 
ral anomaly even at a temperature as low as 10 MeV. It was 
shown that an isotropic and translationally invariant initial 
state of the standard model plasma in thermal equilibrium 
can become unstable in the presence of the global charges[24]. 
The most general form of the polarization operator IT - 7 can 
be written as: 

n 1J (k) = ( k 2 S ij - k i k j '\ IIi(fc 2 ) + ie ijk k k U 2 (k 2 ) (4) 

where, k is a wave-vector and k 2 = |k| 2 . This equation satisfy 
the transversality condition kilT J =0. It should be noted 
here that the Chern-Simon term is oc Y ■ dY , whereas the 
kinetic term is oc (dY ) 2 and therefore the Chern-Simon term 
can dominate over the large length scales. Thus, a non-zero 
value of II 2 when k —Y 0 implying the presence of Chern-Simon 
term in the expression for the free energy. Using the field 
theoretic framework in [24] it was shown that for a sufficiently 
small k < II 2 (A: 2 )/IIi(fc 2 ), the polarization tensor IT"' has a 
negative eigenvalue and the corresponding eigen mode gives 
the instability. 

Recently there has been an interesting development in in¬ 
corporating the parity-violating effects into a kinetic theory 
formalism (see Ref.[31, 36-38]). In this approach the kinetic 
(Vlasov) equation is modified by including the Berry curva¬ 
ture term which takes into account chirality of the particles. 
The modified kinetic equation is consistent with the anomaly 
equation (1). Incorporation of the parity odd physics in ki¬ 
netic theory leads to a redefinition of the Poisson brackets 
which includes contribution from the Berry connection. The 
confidence that the new kinetic equation captures the proper 
physics stems from the fact that the equation is consistent 
with the anomaly equation ( 1 ) and it also reproduces some 
of the known results obtained using the quantum field the¬ 
ory with the parity odd interaction [39]. In fact the “classi¬ 
cal” kinetic equation can reproduce, in the leading order in 
the hard dense loop approximation, the parity-odd correla¬ 
tion of the underlying quantum field theory [30, 40]. The 
modified kinetic equation can also be derived from the Dirac 
Hamiltonian by performing semiclassical Foldy-Wouthuysen 
diagonalization [41, 42]. The modified kinetic equation can 
be applied to both the high density or high temperature 
regimes[41]. Further in Ref. [39] normal modes of the chi¬ 
ral plasma were analyzed using the modified kinetic theory 
in the context of heavy-ion collisions. In this work authors 
have found that in the collisionless limit transverse branch of 
the dispersion relation can become unstable with typical wave 
number k ~ a y/iv, where a is the coupling constant and y, 
refers to the chiral chemical potential. 

Here we would like to note that the authors in Ref. [15] 
have used a heuristically written kinetic equation which is 
consistent with the equation ( 1 ) to study generation of the 
primordial magnetic field. In addition authors have used ex¬ 
pression for the current by incorporating standard electric 
resistivity and the chiral-magnetic effect. The chiral-vorticity 
effect was not considered. It should be emphasized that the 


forms of the kinetic equation used in Ref. [15] and the Berry 
curvature modified theory[37] are very different. As both the 
approaches describe the same physics, it would be interesting 
to see under what conditions they give the similar predic¬ 
tions. Keeping the above discussion in mind we believe that 
it would be highly useful to consider the problem of genera¬ 
tion of primordial magnetic field in presence of the Abelian 
anomaly by using the Berry curvature modified kinetic the¬ 
ory. In this work we incorporate the effect of collisions in the 
modified kinetic theory and derive expressions for the electric 
and magnetic resistivities. The new kinetic framework also al¬ 
lows us to calculate the generation of the primordial magnetic 
field and vorticity. Further in our calculation, we considered 
an isotropic and homogeneous initial state of the particle dis¬ 
tribution function. The magnetic field is generated by the 
unstable transverse modes in presence of chiral charges (Q 5 ). 
This can be seen by integrating equation (1) over space. One 

can write: do 


Q 5 H- H = 0, where, Q 5 = 

7T 
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the helicity T-L = 


V . 


d/x(Y ■ By). The finite helicity state 


can be created even if the initial state has T-L = 0 but Q$ ^ 0. 
Thus the magnetogenesis by net non-zero chiral-charges may 
not require any pre-existing seed field. 

The manuscript is organized into three sections. In 
section-I we briefly state the (3+1) formalism of Thorne and 
Macadonald [43] and the kinetic theory with the Berry cur¬ 
vature. In section-II we apply this formalism to study the 
primordial magnetic field generation in presence of the chi¬ 
ral asymmetry. We also calculated the vorticity generation in 
the plasma due to the chiral imbalance. Section-Ill, contains 
results and a brief discussion. We have shown that our esti¬ 
mated value of the peak magnetic field actually falls within 
the constraints obtained current observations. 


I. BASIC FORMALISM 
Maxwell’s equations in the expanding universe 


In this work we shall study the generation of primordial 
magnetic field at the time when temperature of the Universe 
was much higher than Tr ~ 80 TeV(much higher than the 
EWPT temperature). We intend to solve the coupled system 
of the modified kinetic and the Maxwell’s equations in the 
expanding universe background. Here we note that we ignore 
the fluctuations in the metric due to the matter perturbation. 
For this one needs to write the underlying equations in a gen¬ 
eral covariant form. Interestingly the techniques developed 
in Refs. [44-47] allows one to write the system of kinetic and 
Maxwell’s equations in the expanding background into the 
form that look similar to their flat space-time form. In this 
formalism the well developed intuition and techniques of the 
flat space-time plasma physics can be exploited to study the 
problem at hand. This can be accomplished by choosing a 
particular set of fiducial observers (FIDO’s) [44] at each point 
of space-time with respect to which all the physical quantities 
including liyper-electric and magnetic fields are measured. A 
line element for the expanding background can be written 
using the Friedmann-Lemaitre-Robertson-Walker metric as 

ds 2 = —dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) 


(5) 
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where x, y and z represent comoving coordinates. Here t is 
the proper time seen by observers at the fixed x, y & z and 
a(t) is scale factor. One can introduce the conformal time y 


using the definition y 


J dt/a 2 (t) to write this metric as: 


only T > Tr regime. Thus we write particle number density 
modified by the Berry term as: 

N = lw)'^ 1 + eB - n *V ( 12 ) 


ds 2 = a 2 (y)(—dy 2 + dx 2 + dy 2 + dz 2 )) (6) 


The (hyper)-electric E p h y , (hyper)-magnetic B p h y and the 
current density J p h y are related to the corresponding fiducial 
quantities by transformations: £ = a 2 E p h y , B = a 2 B p h y , 
O' = a 3 J p h y - One can now write the Maxwell’s equations in 
the fiducial frame as: 


V x £ = 0 


dB 

dy 

V ■ £ = 47T p e 

V • B = 0 

V X B = 47 xj + 


d£ 

dy 


(7) 

( 8 ) 
(9) 

( 10 ) 


where B, £, p e and are respectively magnetic field, electric 
field, charge density and current density seen by the fiducial 
observer. 


Kinetic theory with Berry curvature 


The charge p e and current in the Maxwell’s equations 
(7-10) can be calculated using the Berry curvature modified 
kinetic equation which is also consistent with the quantum 
anomaly equation. The modified kinetic equation is given by 


df_ 1 
dy 1 + eflp ■ B 


(e£ + ev x B + e 2 (£ ■ B)Sl p ) ■ ^ 


+ (v + eE x f l p + e(v ■ £l p )B) 


dj_ 

dr 



coll 

( 11 ) 


where e is charge of the particles and has relation with the 
electroweak(EW) mixing angle 9 W \ e = g'cos9 w . Also the 
g' is related with the U(1)y gauge coupling constant a as 
a = g l2 cos 2 9 w /4iv. v = de p /dp = v and e£ = e£ — dt p /dr. 
£lp = ±p/(2p 3 ) is Berry curvature. e p is defined as e p = p(l — 
eB- n p ) with p = \p\- The positive sign corresponds to right- 
handed fermions where as the negative sign is for left-handed 
ones. In the absence of Berry correction i.e. f l p = 0, above 
equation reduces to the Vlasov equation when the collision 
term on the right hand side of equation (11) is absent. 

As we have already stated that we are interested in tem¬ 
perature regime T > Tr Tew ■ At these temperatures 
electrons are massless. The only processes that can change 
electron chirality is its Yukawa interaction with Higgs But at 
this temperature this interaction is not strong enough to alter 
electron chirality. It is here important to note that for tem¬ 
perature smaller than Tr electron mass plays major role in 
left-right asymmetry. Recently in Ref. [48] it has been shown 
that at temperature of the order of MeV, the mass of the 
electron plays an important role in determining the magnetic 
properties of the proto-neutron star by suppressing the chiral 
charge density during the core collapse of supernova. However 
for the present case we ignore the electron mass by considering 


Above equation (11) can be converted to following form by 
multiplying by (1 + eB ■ Sl p ) and integrating over p 


dN „ ^ 2 

~a -h V ■ JT — —e 
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In equation (13) £[ is total current and is defined as = 
E ai 7 a - Here index a denotes current contribution from dif¬ 
ferent species of the fermion e.g. right-left particles and there 
antiparticles. J a is defined as: 


Ja = 
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d 3 p 


a ^ f cl a / & fa \ a rp 

^ + e 


(2tt) 3 
+ e P n p x 


+ e a {£ x <r a ) 


(14) 


where <r“ = / d 3 p/(2Tv) 3 n p f a and e p = p( 1 — e“B ■ f2 p ) with 
p = |p|. Depending on the species, charge e, energy of the 
particles e p , Berry curvature S"2 P and form of the distribution 
function / changes. For the right-handed particle a = R and 
hyper-charge e and for right-handed anti-particle a = R and 
for charge this case is -e etc. It is clear from equation (13) 
that in presence of external electric and magnetic fields the 
chiral current no longer conserved. The first term in equation 
(14) is usual current equivalent to the kinetic theory and re¬ 
maining second and third terms are the current contribution 
by Berry correction. The last term is due to the anomalous 
Hall effect and it vanishes for a spherically symmetric distri¬ 
bution function. If we follow the power counting scheme as in 
[31] i.e. Yp = O(e) , d r = 0(6) (where Y p represents U(1)y 
gauge field) and considering only terms of the order of 0(e5) 
in equation (11), we will have: 


coll. 

(15) 

where, we have taken v = p/p. In the subsequent discussion 
we shall apply this equation to study the evolution of the 
primordial magnetic field. 


d d\ 
dy +V 'dr) 


fa+[ e“5+e“(«xB)-^ 


dfa 

dp 


9fa 

dy 


Current and polarization tensor for chiral plasma 


Here we assume that the plasma of the standard particles 
is in a state of ‘thermal-equilibrium’ at temperature T > Tr 
and at these temperatures the masses of the plasma particles 
can be ignored. We also assume that there exist a left-right 
asymmetry and there is no large-scale electromagnetic field. 
Thus the equilibrium plasma considered to be in a homoge¬ 
neous and isotropic state which is similar to the assumptions 
made in Ref. [15, 24], For a homogeneous and isotropic con¬ 
ducting plasma in thermal equilibrium, distribution function 
for different species are: 

foa(p) = -- ( 16 ) 

exp( P T ) + 1 
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If Sf R and 5fn are fluctuations in the distribution functions of 
the right electron and right-antiparticles around there equilib¬ 
rium distribution. Then we can write perturbed distribution 
functions as 


fR(r,P,v) = fon(p) + Sfn(r,p,r)) (17) 

f R {r,P,v) = foR(p) + 8f&(r,p,v) (18) 

Subtracting equation for a = R from a = R using equation 
(15) one can write: 


d d \ . 

d^ + v 'd^ 


d{B-fl f 


dr 
9f(r,p, rj) 
dp 


+ (£■<# 

op dp 


(19) 


where, f{r,p,r /) = (f R - f n ) and f 0 = foR + f 0 R- Here we 
d f° df° 

have used —— = v^—. This equation relates the fluctuations 
dp dp 

of the distribution functions of the charged particles with the 
induced gauge field fluctuations. The gauge field fluctuations 
can be seen from the Maxwell’s electromagnetic equations (7- 
10). Under the relaxation time approximation, the collision 
term can be written as ( dfa/dr/) co u . « —v c {fa - foa ) (one 
can also look at some studies in chiral kinetic theory with 
collision in Ref. [49, 50]. Next, we take Fourier transform 
of the all perturbed quantities namely £, B and f{r,p,rf) by 
considering the spatio-temporal variation of these quantities 
as exp[— i(ujr) — k ■ r)]. Then using equation (19) one can get 


fk,uj — 


-e[(v ■ £ k ) + v ■ 13 k ){k ■ v 


1 4£o 
J dp 


i(k ■ v — uj — iv c ) 


( 20 ) 


So current contribution for the right handed particle and 
right handed antiparticles in terms of mode function can be 
written using equation (14) as (ignored anomalous hall cur¬ 
rent part): 


3 k R —e 
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In the similar way we can get current contribution from left 
handed particle and left handed antiparticle. So we can ob¬ 
tain total current by adding the contributions from both left 
and right handed particles and antiparticles by putting per¬ 
turbations fkui for all species in equation (21) and adding as: 


J ^n^'dp^ 0R ~ f° R ~ f° L ~ g2 J 

f Ze. 

J 2 tt= 


dp d 


if or + 


foR + foL + /on)> hi — e 2 


27t 2 dp 
if or ~ foR ~ foL + /on)- 


Expression for the polarization tensor I V° can be ob¬ 
tained from equation (22) by writing total current in the 
following form Jy. = Tl lJ (k)Yj(k) using £k = —iujYk and 
Sfc = i(k x Yfe). One can express IT J in terms of longi¬ 
tudinal Pf} — kfy /fc 2 , transverse P%! = (<5 IJ — k r k?/k 2 ) 
and the axial Pf) = ie^ k k k projection operators as Il lfc = 
II lPl? + n tPt + II aPa ■ After performing the angular in¬ 
tegrations in equation (22) one obtains II l, IIt and IIa as 
given below 


U L =-m 2 D ^-[l-uj’L{k)], 


,-r 2 OJUJ 

n T = m D — 


n a = -Er 


1+ ^i m 

,/2 


, uj , T(J , UJ UJ 
1 - - p-) L ( fc ) - -jfc5- 


(23) 

(24) 

(25) 


where, L{k ) = an< ^ u'—u + Wc- Also 

m 2 D = 47ra' + ilR ^ L j and h 2 D = 2a 7r Af * . We have de¬ 

fined Ap = [pr — pl )• In above integrals, we have replaced, 
e by a using relation e 2 = Ana . First, consider case when 
v c = 0. In the limit uj — 1 0, If l and IIt vanish and the parity 
odd part of the polarization tensor II a = h 2 D / 2 tts a'Ap/n. 
Here it should be noted that IIa does not get thermal cor¬ 
rection. This could be due to the fact that origin of Ha 
term is related with the axial anomaly and it is well known 
that anomaly does not receive any thermal correction [51-53]. 
This form of Ha is similar to the result obtained in [24] using 
quantum field theoretic arguments at T < 40 GeV. But in 
the kinetic theory approach presented here such no assump¬ 
tion is made. Normal modes for the plasma can be obtained 
by using expressions for IIt, IIt and IIa. Using equation 
d v F IMV = —AnJ^, we can write to the following relation 

\M~ 1 ) ij Yj{k) = -4 toJI (26) 

where, [M” 1 ] 8 *’ = [(fc 2 — uj 2 )8 lk — k l k k + n* fe ]. Dispersion 
relations can be obtained from the poles of [M~ 1 ]’ ,k , which 
are as given below 

uj 2 = n L , 

uj 2 = k 2 + n T (k)±kn A . 


3 k — 


2 I d£l 
- m D 


v{v .£ k ) 


, f: 

J 47t i(k.v — uj — iu c ) 

h 2 D f dD , „ u 

— J B k -v{v.Bk )} 

d£l ( v X k){v.Bk){k.v) 


Wd 

4 

„2 


/ 

2 J 4n 


47r ( k.v — uj — ii/ c ) 

v{v.B k ){k.v) — (v X k){v.£ k ) 
{k.v — uj — iu c ) 


+ Bk 


( 22 ) 

Here D represent angular integrals. In equation (22), we have 
defined m 2 D = e 2 J + foR + foL + fot), <i = 


One can study the normal modes of the chiral plasma and 
instabilities using these dispersion relations. However, it is 
more instructive to study dynamical evolution of magnetic 
field by explicitly writing time dependent Maxwell equations. 


II. GENERATION OF THE PRIMORDIAL 
MAGNETIC FIELD AND VORTICITY 

Plasmas with chirality imbalance are known to have in¬ 
stabilities that can generate magnetic fields in two different 
regimes: (i) collision dominated k, uj v c [15] and (ii) col¬ 
lisionless case i.e. v c = 0 [39]. In this section we analyze 
how the magnetic fields evolve in the plasma due to these 
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instabilities, within the modified kinetic theory frame work. 
Expression for the total current described by Eq.(22) can be 
written as J k = a Z g£k J +crgB k where <jg and a X g are electrical 
and magnetic conductivities. The integrals involved in equa¬ 
tion (22) are rather easy to evaluate in the limit k, to <C v c 
and one can write the expression for a// and cr l g as: 


a 


ij 

E 


^ i Of. A]l 
3 V c 3v c 7T 


ij £ a ^f l 

B ~ 3 7T 


(27) 

(28) 


Here, we would like to note that the Berry curvature correc¬ 
tion in the kinetic equation gives us an additional contribu¬ 
tion in the expression for <jg which was not incorporated in 
Ref. [15]. First term is the usual dissipative part of the electric 
current and it contributes to the Joules dissipation. The sec¬ 
ond term is due to the chiral imbalance and it does not give 
any contribution to the Joules heating. As we shall demon¬ 
strate later this term is responsible for the vorticity current 
[33]. One can write the total current as = cr//£ k + cr^B^ 
and the Maxwell’s equation: i(k X 13k) 1 = 4-7T J k - Here we 
have dropped the displacement current term (this is valid 

when -- <C 1). Next by taking vector product of k with the 

4na 

above Maxwell equation one obtains (using Maxwell’s equa¬ 
tions and after some simplification): 


dB k 


3v c 


dr] \ 4-7rm; 


k 2 B k - 


a'A/i 


(k x (k x £ k )) 


+ i^^(fcxB fe ) = 0. (29) 


This is the magnetic diffusivity equation for the chiral 
plasma. By replacing (k X £ k ) by —r k in equation 

(29), we can solve this equation without a loss of general¬ 
ity by considering the propagation vector k in z— direction 
and the magnetic field having components perpendicular to 
2 — axis. After defining two new variables: B k = ( B k + iB 2 ) 
and B' k = ( B k — iB k ) one can rewrite equation (29) as, 


The instability can also occur in the collisionless regime 
(v c = 0) [39]. In the quasi-static limit i.e. to -C k one can 
define the electric conductivity as tr| ks ir(mg/2k)S lJ and 
magnetic conductivity <jg ~ (h 2 D /2)8 l: ’. Here it should be 
noted that the above conductivities do not depend upon the 
collision frequency. Similar to the previous case one can take 
the propagation vector in 2 —direction and consider compo¬ 
nents of the magnetic field in the direction perpendicular 2 — 
axis. One can write a set of decoupled equations describing 
the evolution of magnetic field using the variables B k and B k 
as: 


dB k 

dr\ 

d&k 

dij 


7„2 4 a'A iik 

K 3 

9 

Tvmjj 

2k 

jJ2 _|_ 4a/A/xfc 

Trm'jj 

2k 


Bk = 0 , 

Bk = 0 . 


(33) 

(34) 


Here we note that if one replaces d/dr] by —ito equations (33) 
and (34) gives the same dispersion relation for the instability 
as discussed in Ref. [39]. 


Vorticity generated from chiral imbalance in the 
plasma 


It would be interesting to see if the instabilities arising due 
to chiral-imbalance can lead to vorticity generation in the 
plasma. In order to study vorticity of the plasma, we define 
the average velocity as: 

< v >=~ ! 7^4 v(Sf R - Sfn + Sf L - Sf z ) (35) 
n J (27t) j 

Here we have used the perturbed distribution function in 
the numerator of the above equation which is due to the fact 
that the equilibrium distribution function is assumed to ho¬ 
mogeneous and isotropic and therefore will not contribute to 
vorticity dynamics. The denominator is total number density 
and is defined as (in Ref. [54] page number 63): 


dB k 

dr] 

dB/ 

dr] 


( ijc_\ 

V 4nrn 'i J 


^2 _ ( 4aVeA|i \ 


V 


(1 + 




(k 2 + (k 

\inm d J V nrn D / 


(1- 


a.' A/i,fc \ 


Bk = 0, 


B/ = 0 . 


(30) 


(31) 


Thus the magnetic field vector B can be decomposed into 
these new variables Bk & B' k having definite helicity (or cir¬ 
cular polarization). The effect of Ohmic decay is already there 
into the above equations due to inclusion of collision frequency 
v c . It should be noted here that if a A]ik/ivm 2 D -C 1 equa¬ 
tion (30) is similar to the magnetic field evolution equation 
considered in Ref.[15]. In this limit equation (31) will give a 
purely damping mode. In this case dispersion relation will be 


^ ^particle ^ antiparticle 

=d r d 3 P ( i _ i \ 

f J 0 (2tv) 3 yi + exp(^) 1 + exp(^)J 


(36) 


2 2 

which in the case of chiral plasma gives n = —T (]ir + ]il)- 
We consider k, to <C v c regime, in this case the perturbed 
distribution function say for the right-handed particles can 
be written as: 

Sh,u,R = --[(»■ e h ) + ±(v ■ B k )(k ■ «)]^A (37) 

v c 2p ap 

If we add the contribution for all the particles species and 
their anti-particles one can write the numerator in equation 
(35) as: 


ito 


3v c 2 _ 4a'v c A]i 
47 t m d Trm 2 D 


(32) 



1 

V c 


(T 2 ]i 2 R + p \ \ 
V 3 2tt 2 ) 


£k 



(38) 


In Appendix we have shown that the dispersion relation we 
have found here using kinetic theory matches with the disper¬ 
sion relation obtained in [39]. 


Above we have neglected —r ^ R r ^ r> ^ L in comparison to one, 
27T 2 T 2 

as the value of ]ir/T and ]il/T are very small (O(10 - )). One 
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can write average velocity as follows 


< v k >= 



1 1 
2u c ]LR + 


£k 


(39) 


Now vorticity can be obtained by taking curl of the equation 
(39) and assuming that the chemical potentials and temper¬ 
ature are constant in space and time: 


< u>k >= i 


7T 3 2v c n R + n L 


(fc x £ k ) 


(40) 


One can find contribution of the vorticity to the total current 
from equations (27) and J kul = a^£k 2 + rr z ^B J k . By using 
equation (40) the vorticity current can be written as 


Jo ~ - A)u = & (4i) 

Thus in absences of any chiral-imbalance there is no vorticity 
current. Here we note that our definition agrees with the 
Ref. [55]. In Appendix we demonstrate that our kinetic theory 
is also consistent with second law of thermodynamics. Further 
using equation (40), one can eliminating (fe X ft) in equation 
(29) and obtain: 


d&k 3u c 

dr/ 


k 2 &k + i 


.y/'klWt'Vc , 2 2 w, „ V 

2 - KBr - Mn)(« X Wfe) 


+ i 


4a' u c 


(]ir — ]iL){k x Bfc) = 0. 


(42) 


In this equation, second term is usual diffusivity term, how¬ 
ever third and fourth terms are additional term. Which re¬ 
spectively represent vorticity and chiral magnetic effects on 
the chiral plasma. Therefore equation (29) actually contain 
terms due to vorticity and magnetic effect. The saturated 
state of the instability can be studied by setting d v Bk = 0 
in equation (42). After taking a dot product of the equation 
(42) with fluid velocity v k after setting d v Bk = 0 one can 
obtain, 


~ ' Bfc = 0. (43) 

Here we have defined 5 = a (hr — ]j,l)/T. We can write 
expression for magnetic field, which satisfies above equation 
(43) as: 


B k 


g(k)k x 


16T5 

u k - i v k 


(44) 


Where g(k ) is any general function, which can be determined 
by substituting the above expression for the magnetic field 
into equation (42) in the case of steady state. In the very 
large length scale i.e. k —> 0: 


ff( fc ) 


_3_ p- y? R _ 

32 V a' 3 (mh - Bl) 2 


(45) 


So for a very large length scale k —> 0, magnetic field in the 
steady state is: 


Sfc = —i\ 


-7-\ 

4a' (]ir — fiL) 


(46) 


This equation relates the vorticity generated during the in¬ 
stability with the magnetic field in the steady state. 

However in the collisionless regime (w < ): and u c = 0), 
one can have an instability described by equation (33) with 


typical scales k ~ a Afi and |w| ~ a' 2 T5 [39]. Using the 
expression for electric and magnetic conductivities for modes 
in this regime one can write the magnetic diffusivity equation 
as: 


d& h 

dg 


k 2 

4nai 


T$ 

Bk — i - (fc x Bk) = 0 

7T(Tl 


(47) 


where ai = nmjj/2k. Here it should be noted that unlike 
equation (42), the above equation does not have a vorticity 
term. The last term on the left-hand side arises due to the 
chiral-magnetic effect. In the steady state ( d v B k = 0), one 
can gets V X B = (4T5)B. This equation resembles the 
case of magnetic field in a force free configuration of the con¬ 
ventional plasma where the plasma pressure is assumed to be 
negligible in comparison with the magnetic pressure [56]. But 
for our case no such assumption about the plasma pressure is 
required. 


III. RESULTS AND DISCUSSION 


In the previous sections we have applied the modified ki¬ 
netic theory in the presence of chiral imbalance and obtained 
equations for the magnetic field generation for both the col¬ 
lision dominated and the collisonless regimes. The instability 
can lead to generation of the magnetic field at the cost of the 
chiral imbalance. This can be seen from the anomaly equa¬ 
tion til — riR + 2c/W=constant above T > 80 TeV. Where 

A 1 l,rT 2 


til,r 


6 


and H is the magnetic helicity defined as: 


H = 


vl 


d 3 x(Y ■ By) 


(48) 


One can estimate strength of the generated magnetic field 
as follows: From equations (1-2) one can notice that right- 
handed electron number density ur changes with the Chern- 
Simon number n cs of the hyper-charge field configuration as 

1 q' 2 

AriR = -y 2 Rn cs - Here ncs « T 7 r-^kY 2 and A n R = ]i R T 2 = 

2 I 67 U 

88 

—— <5T 3 [54[. From this, one can estimate magnitude of the 
783 

generated physical magnetic field to be 


d phy 

jd y 


2 k8 


2 a'T 


rj~iZ 


(49) 


where we have used By ly ~ kY and k 1 is physical length 
scale, which is related with the comoving length by k~ ky = 
(a/kc) 1 - 

Now consider the regime w, k <S v c where dynamics for the 
magnetic field is described by Eqs. (30-31). Eq.(30) clearly 

gives unstable modes for ( - tt | k < 1 is satisfied. How- 

\nm 2 D J 

ever, Eq.(31) gives a purely damping mode if the condition 


\nmjj 

tion as 


^ k <S 1 is satisfied. 
\3nav c J \ 3-7T 


One can rewrite this condi- 



k <C l.Here we have used 


m 2 D = 3u c a with u c ~ a' 2 ln{^-)T [57] and a = 100T. Thus 

a' 

for k -C v c and 5 <C 1, equation (31) can only give purely 
damped modes. For these values of k and 5, equation (30) 
assumes similar form as the equation for the magnetic field 
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dynamics considered in Ref. [15]. If one replaces — by —icj 

in Eq.(30), the dispersion relation for the unstable modes can 
be obtained. The fastest growth of the perturbation occurs 
ST8 

for Umax i ~ —— and the maximum growth rate can be found 
16 T 2 5 2 

to be Ti ~-=— v c . Here we note that our kmax l differs 

37t mj, 

by a numerical factor from the value of k where the peak in 
the magnetic energy calculated using chiral magnetohydro¬ 
dynamics [23]. For 8 ~ 10~ 6 and a ~ 10 -2 one can show 

that ^ maxl i anc i —F l is satisfied. For these values of 

Vc Vc 

kmax i, ol and 8 one can estimate magnitude of the generated 
magnetic field using equation (49). We find B ~ 10 26 Gauss 
for a ~ 10~ 2 and the typical length scale A ~ 10 J /T. Here we 
would like to note that typical Hubble length scale ~ 10 13 /T 
which is much larger than the typical length scale of instabil¬ 
ity. Our estimate of magnetic field strength B in the collision 
dominated regime broadly agree with Ref. [15]. Here we note 
again that equation (30) includes effect of the Ohmic decay 
due to presence of the collision term. Our analysis shows 
that Ohmic decay is not important for the instability. Fur¬ 
ther we have shown that the chiral instability can also lead 
to generation of vorticity in the collision dominated regime. 
Typical length scale for vorticity is similar to that of the mag¬ 
netic held. From Eq.(46) magnitude of the vorticity to be 
cox ~ 10 _4 B/T. 

Next, we analyze the chiral instability in collisionless 
regime -C oj k, considering equations (33-34). Here one 

S8T 

finds the wave number kmax 2 = —— at which the maximum 


growth rate T 2 = 


and 


1 T 3 8 3 
2n m 2 n 
3 S 3 


Now 


8_5_ 

9 a' 2 


< 1 


£2_ 

u c 8n 2 a 3 ln(l/a ) 

the allowed values of 8. For 


1 this puts constraints on 


1CT 


10" 2 and 


T ~ Tr, one can estimate magnitude of the magnetic held 
to be 10 31 Gauss. Typical length scale for the magnetic held 
A 2 ~ 10/T and which is much smaller than the length scale in 
the collision dominated case. This is expected as the typical 
length scale associated with kinetic theory are smaller than 
the hydrodynamical case (related with the collision dominated 
regime). 


The upper and lower bounds on the present observed mag¬ 
netic held strength from PLANCK 2015 results [58] and 
blazars [59, 60] are between 10~ 17 G — 10 _9 G. However, re¬ 
cently in Ref. [61] it has been shown that if the magnetic held 
is helical and created before the electroweak phase transi¬ 
tion then it can produce some baryon asymmetry. This can 
put more stringent bounds on the magnetic held (10 _14 G — 
10~ 12 G). Since the magnetic helds and the plasma evolutions 
are coupled, the produced magnetic held may not evolve adia- 
batically i.e. like a(rj)~ 2 due to the plasma processes like tur¬ 
bulence. Similarly the magnetic correlation length As oc km\x 
may not be proportional to a{r/). Typical values of As for the 
collision dominated and and collisionless cases in our case are 
10 J /T and 10/T respectively. The length scale of turbulence 

can be written as At ~ — = 77 ~ — - H -1 , 

y/e ch + p ch ^e ch + p ch 

where B p is the physical value of magnetic held and e ch 
and p ch are respectively energy and pressure densities of the 
charged particles. g^iT) and gl° tal (T) are the number of de¬ 


gree of freedom of the U( 1) charged particles in the thermal 
bath. For As At the effect of turbulence can be negligi¬ 
ble. However, the maximum value of the magnetic Held (for 
i/ c = 0) is about 10 31 G in our case and this gives At ~ 10 s /T. 
Thus we have As <S At and following Ref. [61] we assume that 
the generated magnetic helds will undergo inverse cascade 
soon after their generation. One can relate B p and As that 
undergoing the process of inverse cascade with the present 
day values of magnetic held Bo and the correlation length Ao 
using the following equations [61]: 


S^(T)= ! 9.3x10» G ( t; X_) 7,3 ( t5 |^) 


2/3 


/ \ \ 1/3 
(104) Gb{t) (50) 

2/3 


As C (T) ~ 2.4 x 10 -29 Mpc 


10 2 GeV 


-5/3 


Ao 


Bo 


10 ~ 14 G 


tc) G ^ T) (51) 


Where Q b (T) = (si° tai (T)/106.75) 1/6 (sf (T)/82.75) 1/6 (g, s 
(T)/106.75) 1/3 and Q x (T) = (gt° tal (T)/106.75)~ 1/3 {gf (T) 
/82.75)~ 1//3 (<?* S (T)/106.75) 1 ^ 3 . From these equations one 
can see that for collision dominated case B p ~ 10 26 G 
can be achieved when Bo — 1CF 12 G and Ao — 100 Kpc. 
However in collisionless regime a value of Bo ~ 10~ 14 G and 
Ao ~ IMpc at temperature T = 80 TeV gives the values 
that we have found in our estimates for the peak value of 
the magnetic held. Thus the values of the magnetic held 
and the correlation length scale estimated by us can be 
consistent with the current bounds obtained from CMB 
observation and necessary for current observed baryon 
assymmetry. Since the value of B p and As for the collision 
dominated case are simillar to that given in the Ref. [62], so 
we believe that they are also consistent with BBN constraints. 


In conclusion we have studied generation of the magnetic 
held due to the anomaly in primordial plasma consisting of 
the standard model particles. We have applied the Berry cur¬ 
vature modihed kinetic theory to study this problem. The 
effect of collision in the kinetic equation was incorporated 
using the relaxation time approximation. We hnd that the 
chiral instability can occur in presence of the dissipation in 
both collision dominated and collisionless regimes. We hnd 
that in the collision dominated case the chiral instability can 
produce a magnetic held of order 10 27 Gauss with the typi¬ 
cal length scale 10 5 /T. These results are in broad agreement 
with Ref. [15]. However in this work authors have used heuris¬ 
tic kinetic equation and the collision term was not explicitly 
written in the kinetic equation. However the expression for 
the total current included the Ohm’s law. We have obtained 
expressions for electric and magnetic conductivities using the 
modihed kinetic theory. We hnd that expression for electric 
conductivity in chiral plasma has a non-dissipative term in ad¬ 
dition to the standard Ohmic term. It is shown that this new 
term is related to the vorticity current term found in the chi¬ 
ral magnetohydrodynamics [63] . Further we have also studied 
the chiral instability in the collisionless regime. It is shown 
that in this regime magnetic held of strength 10 31 Gauss can 
be generated at length scale 10/T. These length scales are 
much smaller than the length scale of the magnetic Held in 
the collision dominated regime. Further the obtained values 
















of magnetic-field strength and the length scale are shown to 
be consistent with the recent constraints from CMB data. We 
have also shown that in the collision dominated regime results 
of kinetic theory agrees with the hydrodynamic treatment. 


Appendix 

In Ref. [55] it was shown that a parity-violating hydrody¬ 
namics violate the second law of thermodynamics d a s a > 0, 
where s“is the entropy current density, unless certain con¬ 
straints on the transport coefficients are imposed. Therefore 
our results in the collision dominated regime should be in 
agreement with Ref. [55]. The most general equations for U(l) 
and entropy currents can be written as [55]: 

v a = v"*+H(n,T)u'* + H B ^T)B a , (52) 

s“=s'“+%,rK + fl B ( ft T)r (53) 

where, i/“ = pU a + aS a , and s' a = sU a — — v a with p and 

U a being the charge density and hydrodynamic four-velocity 
respectively. Entropy density can be found using the thermo¬ 
dynamic relation sT + pp = (e + p) where, e denotes energy 
density and p denotes the pressure. In the collision dominated 
(MHD) limit there is no charge separation in the plasma and 
one can regard total charge density as zero. Using the energy- 
momentum tensor T al3 = (e + p)m“u /3 + pg a 13 one can write 
equation of motion: d a T al3 = j 1 and divergence of the 
entropy current as 

d a (s a - Doj a -D B B a ) = ~{u a - £u;“ - £ B B a ) 

Z)-C±E a B a . (54) 

According to Ref. [55] the second law of thermodynamic is 
satisfied if the following four equations are satisfied: 


one can write F a13 = SF al3 . Also the energy density e, pres¬ 
sure p and flow velocity U a can be written in this scheme as: 
e = £o + Se, p = po + 8p and U a = Uq + Su a respectively. 
Where all quantities with subscript “0” are background val¬ 
ues. The background and perturbed velocities are defined 
as Uq = (1,0, 0,0) and 8u a = (0,<5u) respectively. Here u 
is the three flow velocity. We assumed that background is 
homogeneous and isotropic. The equation of motion for the 
background gives: do£o = 0 & dopo = 0. Since the £o and 
Po are functions of temperature and chemical-potential, we 
regard background temperature as constant. In this scheme 
the zeroth and i-th components of the equation of motion can 
be written as: 

do Se + (eo + po)diSu l 2 3 = 0 (59) 

(eo + po)doSu + g lj dj5p = 0. (60) 

Ignoring the time derivative term in the Maxwell’s equation 
one can write V x SB = 4nSj and by using using expression for 
the perturbed current one can obtain the following dispersion 
relation 

iw = — ± (61) 

47TCT0 CTO 

where we have used k = kz and <5Bfc = 5Bk x x + 5Bk y y. 
It should be noted that Eq.(61) matches with the dispersion 
relation obtained by kinetic theory approach [equation (32)]. 

We first emphasize that there is no current in the back¬ 
ground and therefore the transport coefficients that appear in 
the expression for the perturbed current depends only on the 
background temperature and the chemical potentials. Now 
consider Eqs. (55-58) which for the background quantities can 
be described by the following two equations: 

-2D bo + t£-=0 (62) 

1 0 

jr-~ C Y 0 = 0 ( 63 ) 


o a D - D - £d a — = 0 

e + p T 

(55) 

d a D B - D b - C B d a ^ = 0 

e+p 1 

(56) 

-2 D b + |=0 

(57) 

— — C— = 0 

T T 

(58) 


In Ref. [55] these equations are solved and one can know 
dependence of £, D and D B on p & T upto an arbitrary 
constant. Next, we assume the perturbation scheme consid¬ 
ered for the kinetic approach. For no background field tensor, 


These equations are satisfied by each species considered. So 
one can write £ 0 = £ R0 + £lo and £bo = £bho + Using 

expression for £o from the kinetic equation, one can calculate 
D b o from Eq.(62) and it agrees with the expression obtained 
in Ref. [55] and also the expression obtained for £b using 
kinetic theory is in agreement with it. Thus we have shown 
that the modified kinetic theory respects the constraint 
implied by the second law of thermodynamics. 
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